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ABSTRACT 



INTEGRAL observations provide a large amount of data 
on accreting binary systems. 

To study the spectral emission of these sources it is neces- 
sary to perform timing analysis and phase resolved spec- 
troscopy, but these tasks are really cumbersome if per- 
formed with the tools based on the imaging extraction 
methods usually used for coded mask instruments. 

Here, we present the software we have developed for the 
ISGRI instrument which allows to extract in a fast way 
light curves, pulse profiles, and phase resolved spectra, 
making data reduction a much easier task. 
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1. INTRODUCTION 



ent energy ranges. Shadowgram deconvolution is how- 
ever a time-consuming process and repeating it on many 
energy bins makes the whole spectral extraction process 
very slow. A much faster way to extract any scientific 
product (light curves, spectra, pulse profiles, etc.) when 
the source position is known, is by using the PIF infor- 
mation. 

For each pixel, the PIF is the fraction of the pixel area il- 
luminated by a source (i.e. PIF=1 if fully illuminated and 
PIF=0 if not illuminated at all). Selected two illumination 
thresholds, PIF on and PIF Q ff , we can define two classes of 
detector pixels, based on their illumination by the source 
of interest: 



ON : PIFy > PIF on 
OFF : PIFy < PIF off 
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Pulse-phase resolved analysis is a powerful method to 
understand the physical emission process in isolated 
or binary pulsars. To perform phase resolved spec- 
troscopy with ISGRI 12] using standard tools it is how- 
ever necessary to follow a cumbersome and considerable 
time-consuming procedure; we have then developed a 
dedicated software that allows to make the data reduc- 
tion in a much faster and easier way. The software uses 
the pixel Photo Illumination Fraction (PIF) to extract sci- 
entific data. 

Spectral extraction methods which use the PIF informa- 
tion have already been implemented for other coded mask 
instruments (e.g. in the official OSA software for JEM-X 
and in the standard extraction software for the BAT in- 
strument on board SWIFT), but to apply this method to 
the ISGRI detector it is necessary to introduce a correc- 
tion step to take into account the effects of its non-ideal 
imaging Point Spread Function (PSF). 



For these two pixel classes, scientific products (light 
curves, spectra, pulse profiles, etc.) can be easily accu- 
mulated, like for a conventional collimated instrument; 
the background subtraction is obtained by the difference 
between the ON and OFF products, after proper normal- 
ization for the number of pixels in each class and, if nec- 
essary, for their efficiency. 

The two illumination thresholds, PIF on and PIF Q ff , can be 
optimized so that the final product will have the best S/N 
ratio. It can be easily verified that for ISGRI, the optimum 
S/N ratio for background dominated sources is obtained 
when PIF on ~ 0.75 and PIF off « 0.25. 

Since in this extraction method the shadowgram decon- 
volution and image analysis steps are avoided, the whole 
procedure is very fast and the processing times are lim- 
ited only by disk I/O operations. 



2. PIF EXTRACTION METHOD 



Spectral extraction in coded mask instruments can be per- 
formed by analysing the sky images generated in differ- 



The outlined extraction method is however strictly valid 
when one single source is present in the instrument Field 
Of View (FOV) or when the imaging Point Spread Func- 
tion (PSF) is ideal. In order to apply this method to 
non-ideal imaging detectors a correction step must be 
performed. 
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Figure 1. ISGRI shadowgram for an almost on-axis 
source. The gaps (two pixels wide) between the eight de- 
tector modules, are evident. 



3. CROSS-CONTAMINATION 

The basic pattern of the IBIS coded mask [3] has been de- 
signed to produce an imaging PSF free of secondary lobes 
for sources observed in the Fully Coded Field Of View 
(FCFOV). The ISGRI detection plane however does not 
fully sample the mask shadow because of the gaps (two 
pixels wide) between the eight detector modules (Fig. [T]>- 

Due to the missing rows and column of pixels in the 
ISGRI detector plane, even when a source is in the 
FCFOV, its PSF has numerous side lobes, forming a cross 
shaped pattern, that can reach intensity levels of several 
percentage of the main peak (Fig. |2). Dead pixels and a 
non-uniform distribution of the pixel efficiencies further 
contribute to the not ideal imaging characteristic of the 
instrument but, if randomly distributed, do not add other 
sidelobes to the PSF. 

A non-ideal PSF implies that scientific products extracted 
by using the simple 'ON/OFF' background subtraction 
method may be affected by the presence of other sources 
in the same FOV It is then necessary to compute and re- 
move the effects of this 'cross-contamination'. 
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Figure 2. ISGRI imaging PSF for an on-axis source. The 
numerous sidelobes forming a cross pattern are gener- 
ated by the missing rows and column of pixel in the de- 
tection plane. The image show only the FCFOV so the 
eight main source ghosts are not visible. 



the sources and the off-diagonal elements, Cjj, the cross- 
contamination between two distinct sources. 

The 'cross-contamination' matrix of an ideal coded mask 
instrument is of course diagonal since the illumination 
patterns of different sources are perfectly uncorrected. 
For ISGRI the off-diagonal elements assume non-zero 
values which depend on the relative source positions; in 
most cases they are of the order of 0.1% but they can 
reach 10% when a source is on one of the main sidelobes 
of the PSF of another source. Higher values can be as- 
sumed when the PSF of a source partially overlaps with 
the PSF of another one (or with its ghost image). 

In any case, the 'cross-contamination' matrix is not sin- 
gular since on each row the highest values are on the 
diagonal (cj^ rj 0.9). The 'decontamination' matrix, 
D = C _1 , allows to compute the 'decontaminated' prod- 
ucts as linear combination of the products generated by 
the ON/OFF background subtraction method. 



4. SOFTWARE VERIFICATION 



3.1. Decontamination procedure 

If the positions of all sources in the instrument FOV is 
known, from their illumination patterns it is possible to 
compute the cross-contamination between any possible 
two of them. This information can be organized in a 
square matrix C (N x N) where N is the number of 
sources in the FOV. The diagonal elements, c,^, of this 
'cross-contamination' matrix are the coding efficiency of 



We have verified that the ISGRI spectra extracted with 
the standard OSA pipeline and with the software we have 
developed, based on the PIF method illustrated above, are 
equivalent in both shape and S/N ratio. In Fig. [3] and 
Fig. |4]the comparison between the Crab spectra and the 
relative S/N obtained with the two methods are shown. 

We have also verified the decontamination procedure by 
extracting the spectrum of a 'fake' source placed in cor- 
respondence of one of the main sidelobes of the PSF of a 



data and folded model 




Figure 3. Comparison between the ISGRI Crab spectra 
extracted with the standard OSA pipeline and with the 
PIF based method. The two spectra are not distinguish- 
able from each other. 
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Figure 5. Crab spectrum extracted with the PIF based 
method (the highest histogram) and spectra extracted in 
correspondence of the most negative sidelobe of its PSF 
before ( the histogram with negative values) and after ( the 
histogram consistent with zero) the application of the de- 
contamination procedure. 



real source. As evident in Fig. [3] where a Crab observa- 
tion has been used, the spectrum of the 'fake' source ob- 
tained without applying the decontamination procedure 
is not consistent with zero; this of course implies that 
the spectrum of a real source in that position would be 
severely altered. 

After applying the decontamination procedure, the spec- 
trum of the 'fake' source becomes consistent with zero, 
thus proving that the cross-contamination problem is 
solved and that the background has been correctly sub- 
tracted. 
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Figure 4. Comparison between the S/N of ISGRI Crab 
spectra extracted with the standard OSA pipeline and 
with the PIF based method illustrated in the text. 



5. PHASE RESOLVED SPECTROSCOPY 



When the source of interest is much stronger than all the 
other sources in the FOV, the cross-contamination can be 
neglected and phase resolved spectroscopy could be per- 
formed using standard software tools (developed for col- 
limated telescopes) on the two events list files correspon- 
dent to the ON and OFF pixel classes. 

In order to make the overall data reduction a much easier 
task than using standard tools, and to apply also the de- 
contamination procedure, we have included in our soft- 
ware also the steps needed to perform phase resolved 
spectroscopy. From the spin period, its derivatives, 
and the eventual binary-orbital ephemeris of a pulsating 
source, the software computes for each event a 'phase 
tag' which is then used to accumulate a two-dimensional 
histogram: the 'phase-energy' matrix (like the one shown 
in Fig- |6j. 

Following the same procedure used to extract spectra and 
light curves, two 'phase-energy' matrices (for the ON and 
OFF pixel classes) are accumulated and then subtracted 
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imaging analysis, the accumulation time, for a typical IN- 
TEGRAL science window, is reduced to the order of few 
seconds (limited by the data I/O operations). 

Moreover, since the processing time is independent on 
the number of energy or time bins used, spectra and light 
curves can be computed at the maximum resolution and 
binned successively, depending of the desired S/N ratio. 
In this way it is not necessary to repeat the whole accumu- 
lation process again (as needed when using the standard 
analysis software) if one should decide to change the size 
of the energy or time bins. 



Figure 6. Example of 'phase-energy' matrix for 
1A 0535+262 observed with ISGRI. The evolution of 
pulse profile with energy is evident. 

after proper normalization; if multiple sources are in the 
same field of view, the decontamination procedure is also 
applied. 

Once the background subtracted 'phase-energy' matrices 
have been generated, phase resolved spectra and pulse 
profiles can be obtained by performing simple selections 
on its rows or columns. In this way, when it becomes 
necessary to change phase selections or energy binning, 
as it can often happen during the analysis of a pulsating 
source, it is not necessary to repeat the whole accumula- 
tion process again. 

Since the accumulation process, even when performed 
with fine phase and energy binning, is very fast (« 1 sec. 
per sew.) and the extraction of any number of spectra or 
pulse profiles from the 'phase-energy' matrix requires a 
negligible amount of time, the overall data extraction for 
phase resolved spectroscopy becomes easy and quick. 



Our software includes also tools dedicated to the phase 
resolved spectroscopy of pulsating sources. In particu- 
lar it is possible to fold the data with any spin period and 
derivatives and, if necessary, correct the event arrival time 
to binary system barycenter. The results of the accumu- 
lation process are background subtracted 'phase-energy' 
matrices from which phase resolved spectra and pulse 
profiles can be later extracted, with a simple selection 
software. 

Example of pulsating sources observed with ISGRI and 
whose phase resolved analysis has been performed with 
our softwarePlare 1A 0535+262 H, 1E1 145.1-6141 H, 
and GX l+4[6]. 
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ment that allows to extract spectra, light curves and per- 
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The extraction method is based on the PIF information, 
but it has been modified to take into account the pos- 
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FOV caused by the non-ideal imaging PSF of the detec- 
tor. We have shown, in fact, that in some cases the cross- 
contamination plays an important role and that it can be 
computed and removed performing simple algebraic op- 
erations. 



[1] Ubertini, et al. 2003, A&A, 411, L131, 2003 

[2] Lebrun, E, et al. 2003, A&A, 411, L141 

[3] Goldwurm, A. et al., A&A, 411, L223, 2003 

[4] Caballero I. et al., these proceedings 

[5] Ferrigno C. et al., these proceedings 

[6] Ferrigno, C, Segreto, A., Santangelo, A., Wilms, J., 
Kreykenbohm, I., Denis, M., & Staubert, R. 2007, 
A&A, 462, 995 



We have verified that the scientific products obtained with 
our software are equivalent to the ones produced with the 
ISGRI standard analysis (OSA) in both shape and S/N ra- 
tio, but they are generated in much less time. In fact, since 
the scientific products are generated by performing sim- 
ple algebraic operations directly on the event list, without authors 
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ABSTRACT 



INTEGRAL observations provide a large amount of data 
on accreting binary systems. 

To study the spectral emission of these sources it is neces- 
sary to perform timing analysis and phase resolved spec- 
troscopy, but these tasks are really cumbersome if per- 
formed with the tools based on the imaging extraction 
methods usually used for coded mask instruments. 

Here, we present the software we have developed for the 
ISGRI instrument which allows to extract in a fast way 
light curves, pulse profiles, and phase resolved spectra, 
making data reduction a much easier task. 



ent energy ranges. Shadowgram deconvolution is how- 
ever a time-consuming process and repeating it on many 
energy bins makes the whole spectral extraction process 
very slow. A much faster way to extract any scientific 
product (light curves, spectra, pulse profiles, etc.) when 
the source position is known, is by using the PIF infor- 
mation. 

For each pixel, the PIF is the fraction of the pixel area il- 
luminated by a source (i.e. PIF=1 if fully illuminated and 
PIF=0 if not illuminated at all). Selected two illumination 
thresholds, PIF on and PIF Q ff , we can define two classes of 
detector pixels, based on their illumination by the source 
of interest: 



> 

(N 
m 



Key words: INTEGRAL; ISGRI. 



1. INTRODUCTION 



ON : PIFij > PIF on 
OFF : PIFij < PIFoff 
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Pulse-phase resolved analysis is a powerful method to 
understand the physical emission process in isolated 
or binary pulsars. To perform phase resolved spec- 
troscopy with ISGRI (?) using standard tools it is how- 
ever necessary to follow a cumbersome and considerable 
time-consuming procedure; we have then developed a 
dedicated software that allows to make the data reduc- 
tion in a much faster and easier way. The software uses 
the pixel Photo Illumination Fraction (PIF) to extract sci- 
entific data. 

Spectral extraction methods which use the PIF informa- 
tion have already been implemented for other coded mask 
instruments (e.g. in the official OSA software for JEM-X 
and in the standard extraction software for the BAT in- 
strument on board SWIFT), but to apply this method to 
the ISGRI detector it is necessary to introduce a correc- 
tion step to take into account the effects of its non-ideal 
imaging Point Spread Function (PSF). 



For these two pixel classes, scientific products (light 
curves, spectra, pulse profiles, etc.) can be easily accu- 
mulated, like for a conventional collimated instrument; 
the background subtraction is obtained by the difference 
between the ON and OFF products, after proper normal- 
ization for the number of pixels in each class and, if nec- 
essary, for their efficiency. 

The two illumination thresholds, PIF on and PIF Q ff , can be 
optimized so that the final product will have the best S/N 
ratio. It can be easily verified that for ISGRI, the optimum 
S/N ratio for background dominated sources is obtained 
when PIF on « 0.75 and PIF off « 0.25. 

Since in this extraction method the shadowgram decon- 
volution and image analysis steps are avoided, the whole 
procedure is very fast and the processing times are lim- 
ited only by disk I/O operations. 



2. PIF EXTRACTION METHOD 



Spectral extraction in coded mask instruments can be per- 
formed by analysing the sky images generated in differ- 



The outlined extraction method is however strictly valid 
when one single source is present in the instrument Field 
Of View (FOV) or when the imaging Point Spread Func- 
tion (PSF) is ideal. In order to apply this method to 
non-ideal imaging detectors a correction step must be 
performed. 
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Figure 1. ISGRI shadowgram for an almost on-axis 
source. The gaps (two pixels wide) between the eight de- 
tector modules, are evident. 

3. CROSS-CONTAMINATION 



The basic pattern of the IBIS coded mask (?) has been 
designed to produce an imaging PSF free of secondary 
lobes for sources observed in the Fully Coded Field Of 
View (FCFOV). The ISGRI detection plane however does 
not fully sample the mask shadow because of the gaps 
(two pixels wide) between the eight detector modules 
(Fig. ??). 

Due to the missing rows and column of pixels in the 
ISGRI detector plane, even when a source is in the 
FCFOV, its PSF has numerous side lobes, forming a cross 
shaped pattern, that can reach intensity levels of several 
percentage of the main peak (Fig. ??). Dead pixels and a 
non-uniform distribution of the pixel efficiencies further 
contribute to the not ideal imaging characteristic of the 
instrument but, if randomly distributed, do not add other 
sidelobes to the PSF. 

A non-ideal PSF implies that scientific products extracted 
by using the simple 'ON/OFF' background subtraction 
method may be affected by the presence of other sources 
in the same FOV It is then necessary to compute and re- 
move the effects of this 'cross-contamination'. 
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Figure 2. ISGRI imaging PSF for an on-axis source. The 
numerous sidelobes forming a cross pattern are gener- 
ated by the missing rows and column of pixel in the de- 
tection plane. The image show only the FCFOV so the 
eight main source ghosts are not visible. 

'cross-contamination' matrix are the coding efficiency of 
the sources and the off-diagonal elements, qj, the cross- 
contamination between two distinct sources. 

The 'cross-contamination' matrix of an ideal coded mask 
instrument is of course diagonal since the illumination 
patterns of different sources are perfectly uncorrelated. 
For ISGRI the off-diagonal elements assume non-zero 
values which depend on the relative source positions; in 
most cases they are of the order of 0.1% but they can 
reach 10% when a source is on one of the main sidelobes 
of the PSF of another source. Higher values can be as- 
sumed when the PSF of a source partially overlaps with 
the PSF of another one (or with its ghost image). 

In any case, the 'cross-contamination' matrix is not sin- 
gular since on each row the highest values are on the 
diagonal (c,^ » 0.9). The 'decontamination' matrix, 
D = C _1 , allows to compute the 'decontaminated' prod- 
ucts as linear combination of the products generated by 
the ON/OFF background subtraction method. 



4. SOFTWARE VERIFICATION 



3.1. Decontamination procedure 

If the positions of all sources in the instrument FOV is 
known, from their illumination patterns it is possible to 
compute the cross-contamination between any possible 
two of them. This information can be organized in a 
square matrix C (N x N) where N is the number of 
sources in the FOV. The diagonal elements, a t i, of this 



We have verified that the ISGRI spectra extracted with 
the standard OSA pipeline and with the software we have 
developed, based on the PIF method illustrated above, are 
equivalent in both shape and S/N ratio. In Fig. ?? and 
Fig. ?? the comparison between the Crab spectra and the 
relative S/N obtained with the two methods are shown. 

We have also verified the decontamination procedure by 
extracting the spectrum of a 'fake' source placed in cor- 



